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Particle number fluctuations are studied in relativistic Bose and Fermi gases. The calculations 
are done within both the grand canonical and canonical ensemble. The fluctuations in the canonical 
ensemble are found to be different from those in the grand canonical one. Effects of quantum 
statistics increase in the grand canonical ensemble for large chemical potential. This is, however, 
not the case in the canonical ensemble. In the limit of large charge density a strongest difference 
between the grand canonical and canonical ensemble results is observed. 
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The statistical models have been successfully used to describe the data on hadron multiplicities in relativistic 
*w . nucleus-nucleus (A+A) collisions (see, e.g., Ref. y| and recent review Q). This has stimulated an investigation of the 
properties of these statistical models. In particular, connections between different statistical ensembles for a system 
of relativistic particles have been intensively discussed. In A+A collisions one prefers to use the grand canonical 
ensemble (GCE) because it is the most convenient one from the technical point of view. The canonical ensemble 
(CE) U S S S S 13 or even the microcanonical ensemble (MCE) have been used in order to describe the pp, 
' pp and e + e~ collisions when a small number of secondary particles are produced. At these conditions the statistical 
| systems are far away from the thermodynamic limit, so that the statistical ensembles are not equivalent, and the 
■ exact charge or both energy and charge conservation laws have to be taken into account. The CE suppression effects 
' for particle multiplicities are well known in the statistical approach to hadron production, e.g., the suppression in a 
production of strange hadrons [(| and antibaryons in small systems, i.e., when the total numbers of strange particles 
' or antibaryons are small (smaller than or equal to 1). The different statistical ensembles are not equivalent for small 
, systems. When the system volume increases, V — > oo, the average quantities in the GCE, CE and MCE become equal, 
' i-e., all ensembles are thermodynamically equivalent. 

The situation is different for the statistical fluctuations. The fluctuations in relativistic syst ems are studied in event 
by event analysis of high energy particle and nuclear collisions (see, e.g., Refs. 0; El Hi El an d references therein). 
In the relativistic system of created particles, only the net charge Q = N + — N^ (e.g., electric charge, baryonic number, 
and strangeness) can be fixed. In the statistical equilibrium an average value of the net charge is fixed in the GCE, or 
. ^ ' exact one in the CE, but N + and iV_ numbers fluctuate in both GCE and CE. 

The particle number fluctuations for the relativistic case in the CE were calculated for the first time in Ref. p3 
^ ■ for the Boltzmann ideal gas with net charge equal to zero. These results were then extended for the CE [TEl Hil Il7j 
\ and MCE [I8I.IT9I ] and compared with the corresponding results in the GCE (see also Ref. |2(j). The particle number 
fluctuations have been found to be suppressed in the CE and MCE in a comparison with the GCE. This suppression 
survives in the limit V — > oo, so that the thermodynamical equivalence of all statistical ensembles refers to the average 
quantities, but does not apply to the fluctuations. 

In the present paper we study the particle number fluctuations in relativistic ideal Bose and Fermi gases for non-zero 
values of the net charge density in the GCE and CE. The paper is organized as follows. In Sections II and III we 
calculate the average values for N + and iV_ and discuss the Bose condensation in relativistic gases. These results are 
not new, and we present them in our paper for completeness. In Section IV we consider the N + and iV_ fluctuations 
in the GCE and study the Bose and Fermi effects for particle number densities. In Section V the same calculations 
and study are repeated within the CE. We compare the GCE and CE results and summarize our consideration in 
Section VI. 
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II. AVERAGE PARTICLE NUMBERS 



The relativistic ideal Bose or Fermi gas can be characterized by the occupation numbers n+ and n~ of the one- 
particle states labeled by momenta p for 'positively charged' particles and 'negatively charged' particles, respectively. 
The GCE average values are pl|: 
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where to is the particle mass, T is the system temperature and p is the chemical potential connected with the conserved 
charge Q: 

Q = (N+) g . c .e. - (N-) g . c . e . = ^(n+) 3 ., e . - $3(n-) fl .c. e . • (2) 

p p 

The parameter 7 in Eq. Q is equal to +1 and —1 for Bose and Fermi statistics, respectively (7 = corresponds to 
the Boltzmann approximation). Each level should be further specified by the projection of a particle spin. Thus, each 
p-level splits into g = 2j + 1 sub-levels. It will be assumed that the p-summation includes all these sub-levels too. In 
the thermodynamic limit the system volume V goes to infinity, and the degeneracy factor g enters explicitly when one 
substitutes the summation over discrete levels by the integration, ^2 p ... = gV(2n 2 )^ 1 J °° p 2 dp ... . The particle 
number densities in the GCE are: 

= (N±)g. c .e. _ Y, P ( n p)g.ce. _ g f°° P 2 dp 
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where m* = m/T, p* = p* /T . To be definite we consider p* > in what follows. This corresponds to non- negative 
values of the system charge density pq = p + — p_ > 0. Results for p* < can be obtained from those with p* > by 
exchanging of -ZV+ and N-. In the Boltzmann approximation (7 = 0) one finds: 

„t3 (gT 3 exp(±u*)/7r 2 , to* < 1 
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gT 3 (m*/27r) 3 / 2 exp[-(m* qp A**)] , m* > 1 



where if 2 is a modified Hankel function. For < p* < to* it follows from Eq. J3J: 

g T 3 ^ 7"- 1 

A>± = —5- to* 2 > K 2 (nm*) exp(±np*) . (5) 

n— 1 

Note that for pF ermi the series expansion in Eq. (JSJ is convergent for all The first term, n = 1, in Eq. (JSJl 
corresponds to the Boltzmann approximation and others, n > 1, are the Bose (7 = 1) or Fermi (7 = —1) statistics 
corrections. For any T and p these correction terms lead to p± ose > p^° ltz and p^_ erml < pBoltz £ Qr particle number 
densities. The Bose enhancement and Fermi suppression factors, 

n Bose n Fermi 

rtBose P± TyF&rvai P± /^\ 

U ± = n Boltz ' K ± = n Boltz ' W 

P± P± 

for different values of to* are shown in Fig. 1 as functions of p* . 

One finds that the largest quantum statistics effects at p* = correspond to the massless particles: 

1 r°° T 2 rl-r ^2, i 

R^(p,* = 0, to* - 0) = i / X = i = C(3) * 1.202 , (7) 

2 7 exp(a;) - 1 ^ n 

i£«"»V = m* = 0) - 5 / TTTT = E^" T~ = iC(3)- 0.902, (8) 

2 7 exp(.T) + 1 ^ n J 4 

where £(fc) is a Riemann zeta function (see Appendix A). Note that the values of njj" and Uq contribute to the net 
charge of the system, but for to = they do not influence the system energy. Therefore, the occupation numbers 
n,Q and Uq become arbitrary, and the ideal Bose gas of charge particles with m = has no clear meaning in the 
thermodynamic limit. In what follows the 'massless' Bose gas of charged particles will be understood as the limit 
to* — > at fixed value of p* = 0. 

For to* 3> 1 using the asymptotic of the K 2 function one finds: 

gi to v-^ wJ r -r- *m g 1 m ' 

P± ~ (27r)3 /2 2^^3/2- exp[-n(TO T A* )] = (2 ^ 3/2 7^3/2(7exp[-(m T A 4 )]) (9) 

where X^i^Li z "/ nfc = Lik(z), is a polylogarithm function (see Appendix A). For z = 1 it equals to the Riemann 
zeta function, Lik(l) = C(^) • The series expansion in Eq. J^J for p+ converges rapidly at p* <t^. to*. In this case it 




FIG. 1: The ratios R+ (left) and R- (right) of particle number densities of bosons (7 = 1) and fermions (7 = —1) to those of 
classical particles (7 = 0) are shown as functions of Two upper solid lines show R± ose at m* = 1, 3. Three lower solid 
lines show R^ erm ' a t m* = 0, 1, 3. The vertical dotted lines fi* = 1, 3 demonstrate the restriction fi* < m* in the Bose gas. 
The crosses at the end of the R+° se and R^ ose lines at fi* = 1 and fi* = 3 correspond to the points of the Bose condensation. 
The crosses at fi* — correspond to the limit m* — > given by Eq. ||7J in the Bose gas. The dashed horizontal line on the left 
corresponds to the maximum value for R + oae given by Eq. dTUI ■ 



is enough to add one term n = 2 to the Boltzmann approximation (n — 1) to describe accurately the Bose or Fermi 
effects. The same is valid for p_ at to* ^> 1 for all fi*. 

The condition fjb* < to* is a general requirement in the Bose gas. At fi* — > to* the Bose enhancement factor R Bose 
reaches its maximum value, and the Bose condensation of positively charged particles starts. This maximum value of 
jlBose a £ _ m * j ncreases w ith. m* and reaches its upper limit, 

max[Rl ose {fi* : m*)] = Rj* ose (/i* = to* -> 00) = Li s/3 (1) = C(3/2) -2.612, (10) 

at to* — > 00 (see Fig. 1, left). For negatively charged particles the Bose enhancement factor reaches its minimal value 
at fM* = to*, 

RBo*e ( * =m * ) ^ ! + exp( _ m * ); (U) 

2K 2 {m*) 

and this value goes to 1 (i.e. to its classical Boltzmann limit) from above at m* — > 00 (see Fig. 1, right). 

The requirement fi* < to* is absent for the Fermi gas and for /i>m one finds (see Eq. (|B6|I in Appendix B): 



Fermi 9 ^ 
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-fi* 3 
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while the density for negatively charged particles can be approximated in this limit with the first two terms from the 
right-hand side of Eq. JSJ): n = 1 corresponds to the Boltzmann approximation pB° ltz and n = 2 gives a small 
(negative) Fermi correction. Thus, one finds that R^ erml goes to zero (see Fig. 1, left) and R^ erml goes to 1 from 
below (see Fig. 1, right) at fi* — > oo: 

Fermi fi* 3 exp{-fi*) Ferml K 2 {2m*) 

3 m* l K2\m*) 2A 2 (m*) 



III. BOSE CONDENSATION 



In a standard non-relativistic picture of the Bose condensation the particle number N is a conserved quantity. If 
the system temperature decreases at fixed particle number density p = N/V, the system chemical potential increases 
and reaches its maximal value at T — Tc- Bose condensation starts and a macroscopic part of the system particles 
- known as the Bose condensate - occupies the lowest momentum state at T < Tc- In a relativistic picture the 
conserved quantity is system charge Q = N + — N- . From Eq. (JSJ) one finds for the dimensionless charge density: 

p Q e A e — - p B ^) = ^ Yl * K 2 (nm*) sinh(r^*) . (14) 

q to oto j v ' tt z ^-^ nm* 

y y n=l 
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Bose condensation starts at the point T = Tq when p — ^ max = m . At this point Eq. 1(14^1 is reduced to: 

PQ = ^ E - K * { n / T c) sinh(n/T c ) , (15) 

n=l 

where Tc = Tc/m. Equation l|15(l can be used to write the Bose condensation temperature Tq as the function of the 
conserved charge density pq. The line of Bose condensation Tc = Tc(pq) given by Eq. l|15f) is shown in Fig. 2 (see 
also Ref. [22( and references therein). 




FIG. 2: The phase diagram of the relativistic ideal Bose gas. The solid line shows Bose condensation temperature as a 
function of the conserved charge density. It is given by Eq. 1151 1 where both quantities are expressed in dimensionless form: 
Tc = Tc/m, pQ = PQ/igm 3 ). At Tc <C 1 the line Tc(pq) is given by non-relativistic approximation 116L while at Tc 3> 1 it is 
described by the ultra-relativistic relation 1171 . The points (p~Q,T) under the solid line correspond to the states of the system 
with non-zero values of the Bose condensate. 



At T c < 1 it follows from Eq. HHJ) : 





PQ - 1^1 l^^n = 1^1 C(3/2). (16) 



This corresponds to a non-relativistic limit. The density of negatively charged particles at T ~ Tc behaves as 
/?_ oc exp(— 2/Tc) and can be neglected in a comparison with p + oc (Tp) 3 / 2 . Under these conditions the charge 
number conservation becomes equivalent just to the (positively charged) particle number conservation. One, therefore, 
recovers from Eq. (|16|) the familiar relation, Tc — 3.313/?^ , between the Bose condensation temperature and particle 
number density known in the non-relativistic statistical mechanics |2l( . At Tc S> 1 from Eq. Ijl5(l one finds: 



2Xg^ 1 T c 
=1 



~ ~ 1/2 

This corresponds to the ultra-relativistic limit and leads to a new relation, Tc — 1.732/9q , between the Bose 
condensation temperature and charge number density. 

For fixed value of the conserved charge density pq the chemical potential is constant p* = m* at T < Tc- The 
positively charged particles have to condensate at the lowest quantum level to preserve a constant value of the positive 
charge density in the system. Therefore, one finds at T < T c : 



PQ = 



Bos 



■{T,p*=m*) - p Bose {T,p* = m*) + p c ° nd {T) , (18) 



where the first two terms in Eq. (|18|) are given by Eq. ©, and p c ° nd is the density of positively charged particles 
at the lowest quantum level (Bose condensate). The behavior of p+ ose , p 30 ^ and p c £ nd above and below the Bose 
condensation temperature are shown in Fig. 3. Note that p+ ose and p 3ose are calculated by Eq. J2| with p* < m* at 
T > T c (the value of p* is defined by the equation p Q = pf ose - p Bose ) and with p* = m* at T < T c {p c ^ nd is given 
by Eq. (|18|) at T < Tc and it equals to zero at T > Tc). The Bose condensation is the 3 rd order phase transition with 
a maximum of specific heat at T — Tc- 
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FIG. 3: The solid lines show the ratios p+/p Q , P- / pq 
relativistic ideal Bose gas. The left and right pictures correspond to the values Tc = 0.5, pq ~ 0.06, and Tc = 10, pq ~ 33.3, 
respectively. Both these (pQ,Tc) points belong to the Bose condensation line in Fig. 2 - the first point lies close to the lower 
left corner and the second point lies close to the upper right corner in Fig. 2. The system presented in the left picture can 
be treated within non-relativistic approximation. In this case, pq ~ p+ ose and p BoBe / p+ ose <C 1, so that negatively charged 
particles can be neglected at T < Tc and charge conservation becomes equivalent to particle number conservation. The system 
presented in the right picture demonstrates the Bose condensation in the ultra-relativistic case: both p+ Bose and p Bose arc 
essentially larger than conserved charge density pQ in the vicinity of the Bose condensation temperature T = Tc- 



IV. PARTICLE NUMBER FLUCTUATIONS IN THE GCE 

The GCE fluctuations of the single- mode occupation numbers are equal to plf : 



i ±2 ) 

L p Ig.c.e. 



L p Ig.c.e.) Ig.c.e. 



l p Ig.c.e. 



"p / g.c.e. 
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= v ± 2 
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(19) 



The fluctuations of the macroscopic observables can be written in terms of the microscopic correlator (Art^ An^) 9 . c . e . , 
where a,/3 are + and/or — , which has a simple form, 



Vp 2 5p k S aP 



(20) 



due to the statistical independence of different quantum levels and different charge states in the GCE. The variances 
of the total number of positively and/or negatively charged particles are equal to: 
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The scaled variance w^ c . e . reads: 
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(22) 



where the thermodynamic limit is assumed, and the p-summation is substituted by the integration similar to Eq. @ ■ 
The scaled variances uf.c° se and w^f f" 11 for different values of m* are shown in Fig. 4 as functions of p* . 
It follows from Eq. (H for 7 = 0, 



+ Boltz 
J g.c.e. 
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i.e. the scaled variances for Boltzmann statistics in the GCE are independent of the chemical potential p,* and equal to 



1 for both the positively and negatively charged particles. The Eq. l]22l) leads to the Bose enhancement, 



and the Fermi suppression to 



a Fermi 



.a Bose 
'g.c.e. 
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< 1, of the particle number fluctuations. 
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FIG. 4: The scaled variances w^ c . e . (left) and w g . c . e . (right) given by Eq. £25) for bosons (7 = 1) and fermions (7 = —1) are 
shown as functions of fi* . The two upper solid lines present 0Jg.c.° 3e for m* = 1, 3. The three lower solid lines present u>^^l rml 
for m* = 0, 1, 3. The vertical dotted lines fi* = 1, 3 demonstrate the restriction fi* < m* in the Bose gas. The crosses at the 
end of the lines for uj~^° S£ at (i* — 1 and fi* = 3 correspond to the points of the Bose condensation, ° se diverges at these 
points. The crosses at = correspond to the limit m* — > given by Eq. 1241 in the Bose gas. 



At \x* = the largest Bose and Fermi effects correspond to the massless particles (see Fig. 4): 
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Using Eqs. (|A5IA6|) one finds from Eq. at ^i* < to*: 
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Note that Eq. ^ is valid for 



,- Fe 
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g.c.e 



for all values of /z* > 0. At to* > 1 one finds from Eq. Prj)l: 

Z ^72- exp[-n(TO* T M*)] 



1 + T E J + IW2 exp[-(n+l)(m* T /z*)] x 
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(27) 



The series expansions in Eq. I|27|l converge rapidly for /x* <C to* — > 00. In this case the term with n = 1 is sufficient 
to describe small Bose or Fermi effects: 



c± c . e -1 + 7 2- 3 / 2 exp[-(m* T M*)] 
The same is valid for negatively charged particles at /1* — > 00: 



^ 2 (2to*) , 

1 + 7 ? ^ exp ~ ^ 
2a 2 (m*J 



(28) 



(29) 



The first terms in Eqs. I|28I29(I correspond to the Boltzmann scaled variances Ij23(l . Therefore, for both positively and 
negatively charged particles, the Bose and Fermi corrections approach to zero as 7exp(— to*) at fi* <C to* — > 00. For 
negatively charged particles, these corrections also tend to zero as 7exp(— fj,*) at fi* — > 00. 

The condition fi* < m* is a general requirement in the Bose gas. At fi* — + to* one the scaled variance oj£ c B e ose 
diverges (see Fig. 4, left). This divergence comes from the contributions of the low momentum modes. Introducing a 
dimensionless parameter S satisfying the conditions m* — /i* -C 6 <C to* one finds: 



Jo 



x 2 dx 
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Therefore, it follows, uj+ c B e ose oc (m* — /i*) x / 2 — > oo, as /i* — > m*. On the other hand, the scaled variance for negative 
Bose particles decreases with pf and reaches its minimum at p* = m*. When jj,* = m* — > oo one finds from Eq. 



, — Bose 
Ug.c.e. 



1 + 2~ 3/2 exp(-2m*) , (31) 



so that LOg, c B e ose approaches to 1 from above as p* = m* — ► oo (see Fig. 4, right). 

The requirement p* < m* is absent in the Fermi gas, and for p* — > oo one finds strong Fermi suppression effects 
(see Fig. 4, left) for positively charged particles (see Eq. (1B14|) in Appendix B): 

Q 

, ,+ Fermi /oo^ 
A* 

The scaled variance for negatively charged Fermi particles increases with p*, and from Eq. (|29(1 . 

- Fermi _ 1 _ ^2 (2m*) /_ *x foo^ 

w ff .c e - 1 2K 2 (m*) 6 Pl M > ' W 
so that <jJg C F e erml approaches to 1 from below at p* — > oo (see Fig. 4, right). 

V. PARTICLE NUMBER FLUCTUATIONS IN THE CE 

In the GCE all possible sets of the occupation numbers {n+,n~} contribute to the partition function. Only the 
average value of the conserved charge Q = —n~) is fixed, (Q} g . c .e. — Q, in the GCE, and (Q) g . c .e. is controlled 

by the chemical potential p*. In the CE an exact charge conservation is imposed. This can be formulated as a 
restriction on permitted sets of the occupation numbers {n^, n~}: only those satisfying the relation, 

AQ = ^ ( An+ - An; ) = > (34) 
p 

contribute to the CE partition function. One proves that this restriction does not change the average quantities in the 
thermodynamic limit, if the average charge in the GCE, (Q) g . c .e., equals the charge Q of the CE (of course, T and V 
values are assumed to be the same in the GCE and CE). In particular, 

(JV+) c .e. - (N+) g . c . e . , (N-)ce. = (AMs.c.e. . (35) 

This is what the thermodynamical equivalence of the CE and GCE means as V — ► oo. This statistical equivalence 
does not apply, however, for the fluctuations, measured in terms of w + and oj~ . The formula i|2U|) for the microscopic 
correlator is modified if we impose the restriction of an exact charge conservation in a form of Eq. I|34|) . One finds (see 
the details in Ref. the CE correlator: 

a 2 a 2„@ 

(An«Ang)c. e . - S pk 8 a p «« 2 - Vp q \* . (36) 
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By means of Eq. (|36|l we obtain: 



(37) 



Comparing Eq. i|36|l and Eq. i|2U[l one notices the changes of the microscopic correlator due to an exact charge conser- 
vation. Namely, in the CE the fluctuations of each mode are reduced, and the (anticorrelations) correlations between 
different modes p ^ k with the (same) different charge states a, (3 appear. These two changes of the microscopic 
correlator result in a suppression of the CE scaled variances in comparison with the GCE ones ujg c e (compare 
Eq. 1)3 7|) and Eq. (|22p). i.e. the fluctuations of both N + and AT_ are always smaller in the CE than those in the GCE. 
A nice feature of Eq. I|37|) is the fact that particle number fluctuations in the CE, being different from those in the 
GCE, are presented in terms of p± and v p 2 given by Eqs. © and (|I9|I . respectively, both quantities calculated in 
terms of (n^) 9 . c . e . within the GCE. 

The Eq. (|I9|I leads to w* 2 = {n p ) g . c .e. in the Boltzmann approximation, so that ^2 p v p 2 = Vp Boltz , and from 
Eq. ij!T7|) one finds (see dashed lines in Figs. 5 and 6): 



± Boltz 



c.e. 



= 1 r >I *\ * = \ \ 1 T tanh (M*)] ■ (38) 

exp(/x*) + exp(— p*) 2 



The Eq. (jSHJ demonstrates the CE suppression effects for particle number fluctuations within the Bolzmann approx- 
imation, e.g., the scaled variances oj^ e Boltz and u>~ e Boltz in the CE at zero net charge density are two times smaller, 



uj+ Boltz = uj- e Boltz = 0.5, than those in the GCE, uj+ B ° ltz = u~ B ° ltz = 1. When the net charge density increases 
the io^ Boltz decreases and tends to at \i* — > oo, while the u>~ Boltz increases and tends to 1. The physical reasons 
of this are seen from Eq. (0J which at p,* ^> 1 gives: p+ ~ pg and p_ <C pq. Therefore, at p* ^> 1 an exact charge 
conservation in the CE keeps JV+ close to its average value Q and makes the fluctuations of N + in the CE small. Under 
the same conditions, (N-) c . e , is much smaller than Q, so that the fluctuations of iV_ are not affected by the CE sup- 
pression effects and they have the Poisson form, as the GCE. The difference between io+ Boltz and ui~ e Boltz , and their 
dependence on p* , are both the new features of the CE. The GCE scaled variances in the Boltzmann approximation 
are equal, u+ B ° ltz = uj~ B ° ltz = 1, and they do not depend on the chemical potential. 

The scaled variances uj^ Bose and uj^ erml , given by Eq. I|37l) . for different values of m* are shown in Figs. 5 and 6 
as functions of p* . At p* = it follows that p+ = /?_ and u+ 2 — v~ 2 . From Eq. (|37[1 we find then for the CE scaled 
variances, 



^ c ± e.(i»* = 0) = ^± c . e .(,i* = 0) . (39) 



According to Eq. (|39[) the CE scaled variances at p* = are two times smaller than the corresponding scaled variances 
in the GCE, e.g., for massless Bose and Fermi particles (see Figs. 5 and 6, and compare with Eqs. (|24I25() '): 

ute. ose (v* = 0, m* - 0) = 1 || ~ 0.684, c±f e ™>* = m* = 0) = 1 || ~ 0.456. (40) 

We study now the CE scaled variances at non-zero values of p* . Let us start with u+ e Bose (Fig. 5, left). At p* — > m* 
it has been found that X^p^ 2 ~~ * 00 ( see Eq. JSOl), thus it follows from Eq. (|3"?|) : 



- »') - - ^"(M- - .,•) x ■ (41) 



The first factor on the right hand side of Eq. (gTJ, uj~ c Bose (p* =m*), reaches its maximum, £(2)/£(3) ~ f .368 
at /i* = m* — > (see Fig. 4, right). When pf = m* — > 0, the second factor on the right hand side of Eq. 14111 . 
p Bose (p* = m*) I p Bose (//* = m*), also increases and goes to 1. Therefore, an upper limit for uj+ e Bose is reached at 
fi* = m* — > (see Fig. 5, left): 



max\u>. 



+ Bose {p*,m*)} = uj- e Bose (p* =0, to*-^0) ~ 1.368. (42) 



c.e. 




FIG. 5: The scaled variances uj^. c B e ose , left, and uj gmC B e ose , right, given by Eq. J37J, are shown as functions of p* . The solid lines 
present ujg B ^ se at m* = 0.01, 0.1, 0.3, 0.5, 1, 2, 3. The vertical dotted lines p* = m* demonstrate the restriction p* < m* 
in the Bose gas. The dashed horizontal line presents a value of £(2)/£(3) — 1.368 which is an upper limit for u!^ e Bo3e reached 
at p* = m* — > (see Eqs. 1421451 h The crosses at p* — m* correspond to the points of Bose condensation. The crosses at 
p* = correspond to u)t. e Boae {p* = 0,m* — » given by Eq. l|iu) . The dashed lines correspond to cj+f° Uz , left and iv~ Boltz , 
right, given by Eq. 
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a 

+ «J 

s 




s 



0.8 



0.6 



— ■ — 1 — 1 1 1 1 1 1 — 


— ■ — 1 — 1 — ■ I ■ I 1 ■ — 








/// 




/// 




/// 


Boltz 

\ 


// / m=2 








m=0 


0.1 


1 



FIG. 6: The scaled variances u>+^ erm! (left) and w c .^ erml (right) are presented by the solid lines for m* = and m* = 2. The 
dashed lines correspond to w+f oltz (left) and olu (right) given by Eq. 



At jj,* = to* — ► oo one finds u> g c Bose (n* — m* — > oo) — * 1 (see Fig. 4, right). Therefore, it follows: 



,+ Bose/,,* _ * 



(/!* = TO* — ► OO) 



p® OSe (H* = TO* -> OO) 1 
pBose^* = TO * _> oo) ~~ C(3/2) 



exp(-2^*) ~ 0.383 exp (-2/x*) , (43) 



so that at /i* = to* — > oo the scaled variance ui^ e Bose goes to zero faster than uj+ e Boltz ~ exp(— 2/j,*). The Fig. 5 (left) 
demonstrates that w + e Bose (/i* = m* = 1) is already smaller than uj^ e Boltz {p* = 1) (Bose 'suppression' !). 
For uj- e Bose ([i* = to*) (see Fig. 5, right) one finds from Eq. Q2J: ' 



^ e Sose (M* - m*) = - to*) , 

so that 

max[uj- Bose {pL*,m*)] = wJ. c * oae (>* = 0, to* -> 0) ~ 1.368 
From Eqs. (|29I44|I it follows, 

w- e Boae ~l + 7 2" 3/2 exp(- 2m*) , 



(44) 
(45) 
(46) 



and uj~ e Bose (n* — m* — -> oo) goes to 1 from above (see Fig. 4, right). 

Now let us turn to the behavior of oj+ e Ferml and ui~ e Ferml (Fig. 6, left and right, respectively). The variance 
J2 P v p 2 f° r the Fermi gas increases as /i* 2 , whereas J2 P v p 2 decreases exponentially, exp(— /i*), for large chemical 
potentials, \x* > 1 . Then it follows from Eq. JSJJl, 



+ Fermi ^ — Fermi v ' — 

W c.e. — w g.c.e. x Fermi 

P+ 



1 X 



m* 2 ivT 2 (m*) exp(-^i*) 
/x* 3 /3 



(47) 



for /i — > oo. Therefore, CL ,+ g i;,er ' mj goes to zero like fx* 3 exp(— /z*) as fi* — ► oo. However, uj+ e Boltz ~ exp(— 2/z*), 
and w + e Ferim becomes larger than uj+ e Boltz (Fermi 'enhancement' !) as /i* — > oo (see Fig. 6, left). Finally, using 
Eqs. (Il:-ii33lljl4l) one finds for u~ e Ferml as n* -> oo: 



, — Fermi _ . , — Fermi . . / i 



1 - ( — ) X 2 (m*)exp(- M *) 



if 2 (2TO*) 

1 " 2i^R) ex P("^ 



Therefore, u; c-e 



Fermi 



goes to 1 at /i* — > oo satisfying the inequalities (see Fig. 6, right): 

, — Fermi ^ , ,— Fermi ^ , — Boltz 
Lu \ Lu \ Lu 

c.e. g.c.e. c.e. 



(48) 



(49) 
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VI. SUMMARY 



The scaled variances for the particle number fluctuations have been systematically studied for the Bose and Fermi 
ideal relativistic gases. The calculations have been done in the grand canonical and canonical ensembles. The analysis 
reveals that in the limit of large (positive) chemical potential the quantum effects and effects of the exact charge 
conservation are absent for negatively charged particles, so that lo~ Bose ~ us~ Ferml ~ Boltz ~ \ m both the GCE 
and CE. However, the strongest quantum effects take place in the limit of large chemical potential for the fluctuations 
of number of positively charged particles in the GCE (see Fig. 4): ^ B ° se — > oo as fi* — ► m* and w£ c F e erml 
as (i* — > oo. On the other hand, just in the limit of large chemical potential we have found the strongest effects 
of the exact charge conservation. The scaled variances w + e Sose and w+ e f crnM (See Figs. 5 and 6) at large fx* are 
very different from those in the GCE. The Bose and Fermi effects in the CE are clearly seen at intermediate fj,*, but 
for n* 3> 1 the effects of the exact charge conservation dominate: oj+ e Bose ~ uj+ e Boltz — * at fj,* = m* — ► oo and 
u)~£~ e Ferml ~ uj^ e Boltz — > at ni — * oo . The summary of analytical results for some limiting values of the scaled 
variances for positively and negatively charged particles in the GCE and CE are presented in Table 1. 







(i* — 0, m* — > 


fi* = m* — > 


H* <C m* — » oo 


H* = m* — ► oo 


m* <C [i* — » oo 




+ Boltz 
Ug.c.e 


1 


1 


1 


1 


1 


Grand 


- Boltz 
Wg.c.e 


1 


1 


1 


1 


1 


Canonical 


. ,+ Bose 
Wg.ee 


1.368 


00 


1 


oo 




Ensemble 


— Bose 
Wg.e.e 


1.368 


1.368 


1 


1 






+ Fermi 
Ug.c.e 


0.912 


0.912 


1 


0.791 







— Fermi 
Wg.c.e 


0.912 


0.912 


1 


1 


1 




,+ Boltz 

W C.£ 


0.5 


0.5 


0.5[1 - tanhOu*)] 


0.5[1 - tanhO*)] -" 


0.5[1 - tanh(ju*)] -> 




,— Boltz 
W c .e 


0.5 


0.5 


0.5[1 + tanh( A t*)] 


0.5[1 +tanh(/i*)] -> 1 


0.5[1 +tanh(/x*)] -> 1 


Canonical 


+ Bose 
Wee 


0.684 


1.368 


0.5 [1 - tanh(/i*)] 







Ensemble 


. .— Bose 
W c .e 


0.684 


1.368 


0.5 [1 + tanh(//)] 


1 






+ Fermi 
W c .e 


0.456 


0.456 


0.5 [1 - tanh(/i*)] 










— Fermi 
W c .e 


0.456 


0.456 


0.5 [l + tanh(/x*)] 


1 


1 



TABLE I: The scaled variances u> + and lo for different statistics in the GCE and CE. The values of fj,* > m* are forbidden in 
the Bose gas. 
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APPENDIX A 

An integral representation of the modified Hankel function K 2 has the form [2j| : 

poo 

K 2 (nm*) — nm*~ 2 \ x 2 dx exp ( — n \J m* 2 + x 2 j . (Al) 

The asymptotic behavior of the K 2 function at large and small arguments is the following: 

K 2 {y) ~ 2y- 2 , y « 1 . (A3) 

Using for z > the expansions (7 = +1,-1,0), 

1 00 00 

— H = E^" lex P(" nz )' 7 7~\ 72 = E^"" 1 " exp[-(n + l)z], (A4) 

exp(z) -7 ^ (exp(^) - 7 ) 2 ^ 
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one finds at [i* < m*: 



x dx 



*2 \- 7 



o exp [(Vm* 2 + x 2 =F A**)] — 7 



V ^— -X 2 [(n + 1) m*] exp [±(n + 1) M *] 
^ — ' n + 1 



'0 [ exp [(Vm* 2 + a; 2 =F A 4 )] ~ 7] n=i 
The polylogarithm function (or Jonquiere's function) is denned as [2 



(A5) 
(A6) 



(A7) 



For z = 1 it equals the Riemann zeta function, 



(A8) 



Some special values of the zeta function used in the paper are: 

=s 2.612, C(2) = y - 1-645, C (3) =s 1.202, C (4) = ^ ^ 1.082 



(A9) 



APPENDIX B 

To obtain the asymptotic expansion of p^ ermt at //* ^> m* one calculates making the variable substitutions and 
integrating by parts: 



x dx 



m* 2 e 



de 1 



n exp(Vm* 2 +x 2 -fi*)+l J m * exp(e- fi*) + 1 3 J m »_^ 
1 
3 



<fo [(y + M*) 2 -m* 2 ] 3 / 2 



exp(y) 



[exp(y) + If 



, \ exp(y) 
* [exp( y) + l] 2 



(Bl) 



The function exp(y)[exp(y) + 1]~ 2 has a maximum at y — and decreases exponentially at y — > ±oo. Expanding the 
function /(y) in a Taylor series at jy = and extending the lower limit of the y-integral in Eq. IjBljl to — oo (this adds 
only exponentially small term proportional to exp(— /i*)) one finds an asymptotic expansion at fx* 3> m*: 



dy f(y) 



exp(jy) 



/m--/i- [exp(y) + l] 2 

It follows for /(0) and its derivatives 

2\3/2 



dy 



/(0) + /'(0) y + - f'(0) y 2 + 



exp(y) 



[exp(y) + l] s 



/(0) = ( M * 2 - m* 2 ) 3 / 2 , /'(0) = 3 M * ( M * 2 - m* 2 ) 1/2 , /"(0) = 

so that /( n )(0) oc (/i*) 3 ~™ as /i* — > oo. The remaining y-integrals are equal to |24j : 

exp(y) . 



2/t 



* 2 



??7 



* 2 



???.' 



2)1/2 



Jo = 



[exp(y) + l] 2 
exp(y) 



[exp(y) + l] 2 



- 2nl J2-jir- = 2 "' (l 2 



fc=i 



(B2) 

(B3) 

(B4) 
(B5) 



for even n = 21, and /„ = for odd n = 21 — 1, (/ = 1, 2, 3, . . .). One finally obtains: 



x dx 



o exp (V m* 2 + a; 2 — /!*) + 1 



(/' 



* 2 



2\3/2 



2/1 



1 

3 

* 2 



/(0) /o + 2 AO) I 2 + 



,* 2 



* 2 



(B6) 
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To find uj + Fermi at n* > ro * one n eeds to calculate integral for a = -7 = 1 (note that W - Fermi ~ w - 
in this limit). Similar to Eqs. l|BllB2jl one finds: 



[ exp (Va; 2 + 
2 
3 



1]' 



f(yo) + /'(w>) (2/ - vo) + 2 /"(wo) (2/ - vo) 2 + 



exp(y) 



[exp(y) + l] 3 



where yo = — In 2 is the point of maximum for the function exp (y) [exp (y) + 1] 3 . The y-integrals in Eq. 
equal to |25| : 

exp(y) 1 



A Q = 


poo 

/ dy 

J —OO 


Ax = 


pOO 

/ dy 

J - — OC 


A 2 = 


p OO 

/ dy 

J —OC 



[exp(y) + l] 3 2 ' 
2/ exp(y) _ _ 1 
2 ' 

^2 



[exp(j/) + l] 3 
2/ 2 exp(y) 



[exp(y) + l] 3 6 



One finds: 



x 2 dx 



[exp (Vm* 2 + x 2 - (x*) + l] 
2 



/(yo) ^0 + f'(yo) (Ai - y A ) + - f"(y ) (A 2 - 2y A, 



vl A ) 



1 * 3 * 2 (* 2 m * 2 

+ (t- — 



Thus for Fermi particles at ^t* — + 00 the following expansions are obtained: 



E« 2 ; 



.9 



VT 3 



g.ce. 



g.c.e. 



27T 2 

5 yr 3 



2tt 2 



3 p V 3 

,.*3 ,,*2 



E< 2 - E( 



lp) g.c.e. 



E< 





M* 






1 










m* 


2 \ , 
1 /i 




2 




gVT 3 



EK: 



g.c.e. 



>yr 3 

2tt 2 



.* 2 



,+ 2 

> /g.c.e. 



.* 2 



27T 2 



(B7) 

(|B7|) are 

(B8) 
(B9) 
(BIO) 



(Bll) 

(B12) 
(B13) 
(B14) 
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